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ABSTRACT 

In  this  thesis  several  incompressible  oscillatory  flow  and 
flutter  problems  were  investigated.  First,  a  previously 
developed  unsteady  panel  code  was  modified  so  that  systematic 
comparisons  with  Theodorsen's  classical  theory  could  be 
accomplished.  It  was  found  that  the  panel  code  is  in  excellent 
agreement  with  the  Theodorsen  results.  Second,  the  panel  code 
was  applied  to  the  analysis  of  bending- torsion  flutter.  Again, 
general  agreement  with  Theodorsen's  flutter  predictions  was 
obtained.  In  the  experimental  part  of  the  thesis  two  flow 
visualization  experiments  were  performed.  First,  the  vortical 
flow  patterns  generated  by  an  airfoil  executing  harmonic 
plunge  oscillations  were  visualized.  In  the  second  experiment, 
the  interference  effects  between  a  stationary  airfoil  and  a 
small  vane  executing  plunge  oscillations  were  explored. 
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I.  INTRODUCTION 


A.  GENERAL 

In  this  thesis,  several  numerical  methods  were  used  to 
analyze  the  flow  about  an  airfoil  performing  unsteady  motion 
in  an  inviscid  incompressible  fluid.  First,  the  unsteady 
motion  of  a  single  airfoil  was  analyzed  after  modifying  the 
U2DIIF  code  [ref .2] .  The  primary  purpose  was  to  verify  the 
code  against  the  proven  theory  of  Theodorsen  for  analyzing  the 
phenomenon  of  flutter.  To  accomplish  this  the  U2DIIF  code  was 
modified  to  calculate  aerodynamic  values  over  a  range  of 
reduced  frequencies  and  then  apply  these  values  to  the 
flutter  analysis. 

Next,  the  propulsive  effects  of  a  plunging  airfoil  were 
verified  through  experimental  methods  using  a  low  speed 
plexiglas  wind  tunnel . 

Finally,  an  exploratory  test  was  conducted  in  the 
department's  smoke  tunnel  to  study  the  interaction  between  a 
plunging  airfoil  and  a  stationary  large  airfoil. 

B.  SCOPE 

Chapter  II  contains  the  modification  of  the  single  airfoil 
U2DIIF  code  into  the  code  UPOT.f  and  extensive  verification  of 
this  code  against  results  produced  by  Theodorsen.  Chapter  III 
describes  the  UPOT  code  and  explains  the  modifications  which 
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were  added  to  solve  the  flutter  determinant.  In  chapter  IV  the 
flow  visualization  experiment  is  described  which  was  performed 
to  study  the  vortical  wake  patterns  produced  by  a  plunging 
airfoil .  In  chapter  V  a  second  experiment  is  described  which 
was  performed  to  explore  a  plunging  airfoil's  potential  for 
control  of  flow  separation. 


2 


II.  SINGLE  AIRFOIL  ANALYSIS 


A.  n2DIFF  PANEL  CODE 

1 .  Geometry 

Figure  2.1  shows  a  representation  of  the  system  that 
is  analyzed  using  the  panel  code.  Shown  are  the  values  for  h 
(plunge)  and  a  (AOA) . 

2 .  U2DIFF 

The  U2DIIF  code  was  developed  by  TENG  [ref  .2]  for  the 
study  of  unsteady  inviscid  and  incompressible  flow  over  a 
single  airfoil.  The  code  is  based  on  the  extension  of  the 
panel  method,  developed  by  Hess  &  Smith  [ref. 4]  for  steady 
potential  flow  problems,  to  include  the  unsteady  motion  of  the 
airfoil  that  is  continuously  shedding  vortices  into  the 
trailing  wake.  This  vortex  shedding  process  is  nonlinear  in 
that  the  wake  vortices  influence  the  flow  over  the  airfoil 
which  in  turn  alters  the  vortex  shedding  as  the  airfoil 
proceeds  in  time. 

The  non-linearity  of  the  unsteady  flow  makes  this 
problem  different  from  the  steady  flow  problem  which  requires 
only  simple  Gaussian  elimination.  Teng  developed  a  code  that 
used  an  iterative  type  of  solution.  Typical  program  output 
includes  the  airfoil  pressure  distribution,  force  and  moment 
coefficients,  and  the  trailing  vortex  wake  pattern.  No 
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Figure  2.1  Airfoil  Geometry 
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distance  between  rotation  point  (elastic  axlsX 
and  nidchord,  positive  if  aft  of  nidchord 


attempt  is  made  here  to  reproduce  the  work  of  Teng  or  to 
explore  the  operation  of  the  U2DIIF  code,  but  the  reader  is 
encouraged  to  review  reference  2 . 

B.  PHASE  PROGRAM 

The  phase  program  was  put  through  some  verification  by 
Neace  [ref. 9]  and  modified  slightly  in  order  to  present 
results  for  harmonic  motion.  The  code  PHV3.f  (phaseshift)  was 
written  by  Neace  to  convert  the  time  dependent  output  of  lift 
and  moment  histories  to  harmonic  output  using  an  iterative 
curve  fit  algorithm: 

F(t)  =Amp*Sin{(iiit+^)  (2.1) 

where  Amp  =  amplitude  of  motion,  u  =  frequency,  and  4)  =  phase 
angle  between  motion  and  the  aerodynamic  forces.  One  primary 
output  of  this  program  was  the  values  of  phaseshift  (4>) 
between  the  AOA  and  coefficients  of  lift  (C^)  and  moment  (C^) 
for  the  pitching  airfoil  and  the  phaseshift  between  the 
plunge  value  (h/2b)  and  the  Cl  and  C,,  for  the  plunging  airfoil. 
The  other  output  was  the  amplitude  of  Cl  and  C,,  for  the 
pitching  or  plunging  case. 

C.  MODIFICATION  OF  t72DIIF  AND  PHASE  PROGRAM 

In  an  attempt  to  make  the  above  mentioned  codes  more  "user 
friendly" ,  the  two  codes  were  combined  into  a  single  code 
named  UPOT.f.  The  modification  involved  a  new  input  file 
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called  UPOT.in  which  gives  the  user  several  options  of 
operation.  The  input  file  can  call  for  the  analysis  of  steady 
flow  only,  straight  and  modified  ramp  motion,  pitch 
oscillation,  plunge  oscillation,  and  the  capability  of 
performing  the  oscillation  analysis  over  a  series  of  reduced 
frequencies.  A  sample  input  file  is  shown  in  Figure  (2.2). 

1 .  Output 

Outputs  from  the  code  have  been  limited  to  reduce  the 
amount  of  computer  space  taken  up  by  the  code  operation.  A 
sample  calculation  was  run  using  the  input  from  Figure  (2.2) 
and  the  on-screen  output  is  shown  in  Figure  (2.3) .  Of  course, 
the  user  can  modify  the  output  portions  of  the  code  to 
minimize  output.  The  following  list  describes  the  input/output 
files  and  the  data  they  contain. 

a.  UPOT.IN:  The  input  file  figure  (2.2). 

b.  CL.d:  This  file  contains  the  various  AOA  values  along 
with  its  corresponding  for  each  time  step. 

c.  CM.D:  This  file  contains  the  various  AOA  values  along 
with  its  corresponding  for  each  time  step. 

d.  PHASE. d:  This  file  contains  the  values  of  non 
dimensional  time  (t)  ,  AOA,  C^,  C„,  for  each  time  step. 

e.  FOR015.DAT:  This  file  contains  the  values  of  non 

dimensional  AOA,  curve  fit  for  C^,  curve  fit  for  C„, 
(used  in  the  phase  portion  of  program) . 

f.  CPSS.d:  This  file  contains  the  steady  state  pressure 
coefficient  for  the  mid  point  locations  of  all  the  air 
foil  panels. 

g.  CPUOOB.d:  This  file  contains  the  unsteady  pressure 
coefficient  for  the  mid  point  location  of  all  airfoil 
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panels,  (in  this  case  the  values  are  for  an  AOA  equal 
to  5  degrees) . 

h.  PHZSWP.d:  This  file  contains  the  phase  information  of 
the  reduced  frequency  sweep  portion  for  the  program. 
The  file  contains  the  phase  angle  of  C^,  and  C„,  and 
the  amplitude  of  C^,  C^. 

i.  FLUTTER. IN:  This  file  contains  information  that  can 
be  used  to  solve  the  flutter  determinant.  It  contains 
Kp,  Re,  Cl  Im,  C„  Re,  Im  for  the  pitch  or  plunge 
case. 
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Page  1 


AIRFOIL  TYPt  :  NACA  0012  AIRFOIL 
SLOWER  =  50  ,  NUPPER  -  50 


I  FLAG  *,  LOWER 
0  50 

AIRFOIL  TYPE 
7 


IRAMP 

0 

FREQ 

.68 

IGOST 

0 

I TRANS 
0 

CYCLE 

2 


lOSCIL 

1 

RFQS7P 

O.Cl 

UGUST 

0. 

DELHX 

.00 

STCYCLE 

60 


ALPI  A 
-3.0 
RFQFNL 

.7 

VGOST 

0. 

DELHY  DELI 
.0  .0 
TOL 
0.005 


naot  X  aoa  values  multiplied  by  1C  <inteqer) 
05  10  20  25  39  50 


Comment  s . . . 

IRAMP  0:  n/a  RFREQ  is  based  on  ful.  chore 

1;  Straight  ramp 
2:  Modified  ramp 

lOSCIL  0:  n/a  RFREQ  is  based  on  full  chord 

1:  Sinusoidal  pitch,  motion  starts  at  min  Aoa 

ITRANS  0:  n/a 

1:  Translational  harmonic  oscillation 

ALPr/ALPMAX  Minimum/MAX  AOA  in  degrees  for  IRAMP/ITRANS/IOSCIL 
MAX  does  not  apply  for  ITRAns 

PIVOT  Location  of  Elastic  Axis  as  a  fraction  of  full  chord 
FREQ  Initial  reduced  frequency  for  program 

RFQSTP  Reduced  freq  step  size  for  a  sweep  of  f req. ' s (enter  0.0  if  only  one  calcu 
latlon  is  desired.) 

RFQFNL  Final  freq  for  the  sweep 

OELHX  Translational  amount  in  the  chordwise  direction  (dist/full  chord) 

DELHY  Max  Translational  amount  in  the  vertical  direction  (h/fullchordfb) ) 

DELI  Min  Translational  amount  in  the  vertical  direction  (h/b) 

CYCLE  :  f  of  cycles  for  oscillatory  motions 

*-In  case  of  ramp,  cycle=1.5  denotes  airfoil  is  held 
at  max  aoa  for  the  duration  of  .5  cycle 
-For  steady  state  solution  set  It  to  0 

NTCYCLE:  f  of  time  steps  for  each  cycle 

CYCLE*NTCYCLE  is  limited  to  200  currently, 

TOL  Tolerance  for  convergence  of  the  unsteady  sol ut ion. (recommend  using  not 
less  than  .001) 

NAOT:  t  of  input  aoa  for  cp  output 

-  angles  should  be  in  increasing  order, 

-  for  oscllatory  motions  angles  should  increase 
first,  then  decrease.  Decreasing  angles  are  for 
the  return  cycle.. 


Figure  2.2  UPOT.IN 


Stdin _  Page  1 
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A:Rb'0:i  TYPE  :  NACA  0012  AIRFOIL 
SLOWER  ^  50  ,  NUPPER  =  50 


I FLAG  (C:NACA,  1: INPUT) 

NO.  PANELS  UPPER  SURFACE 

NO.  PANELS  LOWER  SURFACE 

=  0 
=  50 

=  50 

7 

OSCILLATORY  MOTION,  lOSCIL 
INITIAL  ANGLE  OF  ATTACK 
FINAL  ANGLE  OF  ATTACK 

REDUCED  FREQ.  FOR  OSCIL 
REDUCED  FREQ.  STEP 

FINAL  REDUCED  FREQ. 

PIVOT  POINT 

=  -3.C00C 

=  3.0000 

=  C.6dOC 

0.0100 
=  0.7000 

=  0.3700 

TOTAL  1  OF  CYCLES 

1  of  TIME  STEPS  PER  CYCLE 
TOLERANCE  FOR  CONVERGENCE 

»  2.0C0C 

=  60 

=  0.0050 

FREO  SWEEP 
FREO  =  0.680000 

STEADY  FLOW  SOLUTION  AT  ALPHA  *  -3.0000CO 


1 

2 

3 

4 

5 

6 
1 
8 
9 

10 

U 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 


0.999507 

0.997535 

0.993600 

0.987718 

0.979910 

0.970208 

0.958651 

0.945283 

0.930159 

0.913336 

0.894883 

0,874870 

0.853379 

0.830493 

0.806302 

0.780903 

0.754395 

0.726883 

0.698476 

0.669285 

0.639427 

0.609018 

0.578179 

0.547031 

0.515698 

0.484302 

0.452969 

0.421821 

0.390982 

0.360573 

0.330715 

0.301524 

0.273117 

0.245605 

0.219097 

0.193698 

0.169507 

0.146621 


-1.116717 

-1.108171 

-1.100064 

-1.092523 

-1.085444 

-1.078820 

-1.072610 

-1.066794 

-1.061323 

-1.056157 

-1.051250 

-1,046545 

-1,041994 

-1.037547 

-1.033154 

-1.028773 

-1,024362 

-1,019887 

-1.015315 

-1.010628 

-1,005803 

-1.000836 
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-0,985017 

-0,979449 

-0.973745 

-0.967910 

-0,961956 

-0.955883 

-0.949692 

-0.943373 

-0.936907 

-0.930250 

-0.923337 

-0.916066 

-0.908273 

-0.899722 


-0.086820 

-0.086820 

-0.086820 

-0.086820 

-0.086820 

-0.086820 

-0.086820 

-0.086820 

-0.086620 

-0.086820 

-0.086820 

-0.086820 

-0.086820 

-0.006820 

-0.086820 

-0.086820 

-0.086820 

-0.086820 

-0.086820 

-0.086820 

-0.086820 

-0.086820 

-0.086820 

-0.086820 

-0.086620 

-0.086820 

-0.086820 

-0.086820 

-0.086620 

-0.086820 

-0.086820 

-0.086820 

-0.066820 

-0,086820 

-0.086820 

-0.086820 

-0.086820 

-0,086820 


0.280753 

0.213074 

0.168200 

0.134381 

0.106525 

0.082353 

0.060678 

0.040815 

0.022307 

0.004851 

-0.011780 

-0.027772 

-0.043274 

-0.058416 

-0.073319 

-0.086096 

-0.102850 

-0.117685 

-O.X32698 

-0.147973 

-0.163603 

-0.179660 

-0.196223 

-0.213352 

-0.231108 

-0.249545 

-0.268710 

-0.288653 

-0.309418 

-0.331061 

-0.353649 

-0.377267 

-0.402032 

-0.428106 

-0.455707 

-0.485152 

-0.516866 

-0.551434 


-0.64B0B4 
-0.887088 
-0.912031 
-0.930386 
-0.945238 
-0.957939 
-0.969166 
-0.979380 
-0.968783 
-0.997572 
-1.005873 
-1.013791 
-1.021408 
-1.028793 
-1.036011 
-1.043118 
-1.050167 
-1  .057206 
-1  .064283 
-1,071435 
-1.078704 
-1,086122 
-1.093720 
-1.101523 
-1.109553 
-1.117830 
-1.126370 
-1.135188 
-1.144298 
-1.153716 
-1.163464 
-1,173570 
-1.184074 
-1,195034 
-1.206527 
-1.216666 
-1.231611 
-1.245566 


Figure  2.3a  UPOT  output 
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-1.260822 

40 
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-1.277774 

41 
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42 
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-1.319145 

43 

0.054717 
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-1.345363 
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-1.311093 

45 
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55 

0.020090 
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57 

0.041349 
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0.219097 
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71 
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32 
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73 
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14 
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75 
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79 
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80 
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82 
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83 

0,726883 
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84 

0.754395 

0.891196  -0.086820  0.014851 

0.992547 

85 
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86 

0,806302 

0.893264  -0.086820  0.025328 

0.981255 

81 

0.830493 

0.894605  -0.086820  0.031258 
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88 

0.853319 

0.896208  -0.086820  0.037145 
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89 
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90 
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0.900461  -0.086820  0.052115 

0.973255 

91 

0,913336 
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92 

0.930159 
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0.963630 

93 

0,945283 

0.910146  -0.086820  0.082565 

0.957828 

94 
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0.915603  -0.086820  0.095332 

0.951141 

95 

0.970208 

0.921351  -0.086820  0.110175 

0.943306 

96 
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0.928123  -0.086820  0.121830 
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97 

0.987718 

0.936086  -0.086820  0.149558 
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98 
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99 

0.997535 

0,957020  -0.086820  0.217400 

0.884647 

100 

0.999507 

0.911841  -0.086820  0.280753 

0.B480B4 

•  •  • 

BEGIN  UNSTEADY  FLOW  SOLUTION 

istep  alpha 

time  r\ltr  cl,  cd, 

,  cm 

1 

-3.0000 

0,0000  I  -0.3479 

0.0002  -0.0403 

2 

-2.9836 

0.1540  0  -0.3372 

•0.0005 

-0.0424 

3 

-2.9344 

0.3080  0  -0.3247 

•0.0012 

■0.0434 

A 

-2.8532 

0.4620  0  -0.3101 

•0.0019 

'0.0441 

Figure  2.3b  UPOT  output 
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0.1653 

0.0024 

0.0280 

35 

2.7406 

5.2360 

‘ 

C.1538 

0.0016 

0.0289 

36 

2-5981 

5.3900 

0 
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Figure  2.3d  UPOT  output 
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Figure  2.3e  UPOT  output 
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Figure  2.3f  TJPOT  output 
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-0.0002 

0.0175 

107 

-0.3136 

16,0874 

0 

-0.0782 

0.0005 

0.0145 

108 

-0.6237 

16.2391 

0 

-O.lOU 

0,0012 

0.0113 

109 

-0.9270 

16.3909 

0 

-0.1230 

0,0020 

0.0079 

110 

-1.2202 

16.5427 

0 

-C,1437 

0.0028 

0.0044 

111 

-1 . 5000 

16.6945 

1 

-0.1630 

0.0035 

0.0009 

112 

-1.7634 

16.6462 

0 

-0.1807 

0,0042 

-0.0026 

113 

-2.0074 

16.9980 

0 

-0.1965 

0.0048 

-0.0061 

114 

-2.2294 

17.1498 

0 

-0,2103 

0.0052 

-0.O096 

115 

-2 .4270 

17.3015 

0 

-0.2219 

0,0055 

-C.0130 

116 

-2.5981 

17.4533 

1 

-0.2312 

0.0056 

-0.0163 

117 

-2.7406 

17.6051 

0 

-0.2382 

0.0055 

-0.0194 

118 

-2.8532 

17.7568 

c 

-0.2427 

0.0053 

-O.C  23 

119 

-2.9344 

17.9086 

0 

-0.2446 

0.0049 

-0.025C 

120 

-2 . 9836 

18.0604 

c 

-0.2441 

0.0044 

-0.0274 

PHASE  SHIFT  ANALTSIS 
FREO  =  0.6900000 


Figure  2.3g  UPOT  output 


stdin 


Pages 


AyPLITunE;  clanp,  cmamp 
PHASE;  Clp,  cmp 

AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  WBAR 


0.J300964 

193,5791 

1.5937672E-03 

-0.2118947 


3.4413978E-02 

-37.85255 

9,4U9796£-02 

-7.2855391E-C3 


Figure  2.3h  UPOT  output 
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2.  Difference  Between  UPOT  and  U2DIIF/Phase 


The  input  file  format  was  changed  along  with  the 
following: 


•  The  program  can  now  analyze  a  pitch,  plunge  or  ramp  motion 
that  starts  from  any  minimiun  value  of  Alpha  or  plunge 
(h/2b)  .  Previously,  the  program  only  accepted  the  initial 
position  of  zero.  This  program  does  not  need  to  go 
through  the  origin. 

•  The  phase  portion  of  the  program  was  changed  to  curve  fit 
Cl  and  Cjj  to  a  cosine  function: 

F(t)  =Ajmp*Cos{(ot+^)  (2.2) 

where.  Amp  =  amplitude  of  motion,  co  =  frequency  of  motion, 
0  =  phase  angle  between  motion  and  the  aerodynamic 
forces.  This  was  done  since  the  alpha  and  plunge  values 
were  allowed  to  start  from  a  new  zero  position. 

•  The  phase  portion  uses  the  middle  180  degrees  of  the  final 
360  degree  cycle  specified  in  the  UPOT. in  file.  This 
change  was  done  to  capture  an  all  positive  area  of  the 
cosine  curve  for  phasing  analysis.  The  program  integrates 
this  portion  of  the  cosine  curve,  and  for  proper  code 
operation  the  area  under  the  curve  must  be  kept  to  one 
sign.  If  the  areas  of  integration  were  chosen  to  include 
both  sides  of  the  axis,  then  the  code  would  produce  errors 
near  90  and  270  degrees. 

3.  UPOT  Verification 

The  code  UPOT  did  not  incorporate  any  drastic  changes 
to  the  prior  codes,  but  the  original  code  had  never  been 
extensively  compared  to  prior  theories  over  a  wide  range  of 
reduced  frequencies.  When  conducting  these  comparisons,  it  is 
easy  to  become  confused.  This  section  will  go  through  the 
comparisons  slowly  to  help  alleviate  that  problem. 


17 


a. 


^panel  •  ^Tbeodorsaa  * 

The  equation  for  reduced  frequency  is: 


,  o)2jb  j,  _  coJb 
U  •=  U 


(2.3) 


where:  w  =  frequency  of  oscillation  (rad/sec) 
b  =!  half  chord  (units  to  match  U) 

U  =  free  stream  velocity  (units  to  match  b) . 

The  difference  between  Kp  and  lies  in  the  fact  that  Kp 

calls  for  the  full  chord  and  calls  for  the  half  chord, 

hence,  it  is  important  to  remember  that  Kp  is  twice 
b.  Aerodynamic  Forces 

The  aerodynamic  forces  problem  of  simple  harmonic 
motion  about  an  equilibrium  position  was  solved  theoretically 
by  Theodorsen  in  NACA  TR-496  [ref. 10]  and  outlined  by  Fung  in 
[ref .5] .  The  complex  equations  were  simplified  using  the 
simple  harmonic  motion  equation  and  resulted  in  the  following: 

L=npb3w2  (L^l  +  [L.-  ( A+a)  LJ  a+  [Lp-  (c-e)  Lj  p)  (2.4) 

M=Tipb*w^(  [M^-(-|+a)  Lj  + 

[M.-(4+a)  (A.+W^)  +(-|+a)2L^]a  (2.5) 

+  [Afp-{-|+a)Lp-(c-e)M2--(c-e)  ( -| +a)  Lj  P) 

L,  M  are  the  lift  and  moment  per  unit  span  of  the 
airfoil  about  the  elastic  axis,  b,  h/b,  a  and  cy  (radians)  ,  are 
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shown  in  Figure  (2.1).,  Lj^,  Lp,  and  are  defined  by 

Scanlan  [ref .6, pp. 412 -424]  for  various  values  of  and  e. 
This  analysis  will  not  cover  airfoil  aileron  combinations. 
Therefore  |S  becomes  zero  and  equations  2.4  and  2.5  reduce  to: 

L=npb^w2  +  [L.  -  ( ^  +a)  LJ  o )  (2.6) 

JD  ^ 


W=npb"<o2(  [/4-(l+a)Lj  ^  + 
[W.-(|+a)  (L,+i%,)  +  (|+a)2Lja)e''“^**^> 


(2.7) 


The  UPOT  panel  code  used  the  following  equations 
in  defining  lift  and  moment: 

Cr  =  =jRl+lle^ (® f+i,)  X  (2.8) 

^  2qb  y  ^ 

C„=  — ^  =Jr^ + j 2 e i (« t**»)  ^  ^tan'i-^  (2.9) 

"  4qi)2  V  ^ 

where  and  are  the  real  and  imaginary  parts  of  C^/  and  Rm 
and  are  the  real  and  imaginary  parts  of  Cjj.  For  the  same 
conditions  the  lift  (L)  and  moment  (M)  should  be  the  same  for 
both  the  panel  code  and  Theodorsen.  This  fact  allows  for 
comparison  of  the  magnitude,  real,  imaginary  and  phase  of  lift 
and  moment. 

For  lift:  Lj.  (eqn  2.6)  equals  Lp  (eqn  2.8) 


L9 


npjb^io^  [L^h/Jb+  (L.-  (Vz+a)  L*)  a]  =2g2yi?|+j|e^‘"‘*^^’  (2 . 10) 


After  canceling  ^ 


npjb^w^  [L^(i3/jb)  +  (L,-  (Vz+a)  L^)  a]  =2qb^Rl+ll 


(2.11) 


For  pitch  case,  h/b  =  0,  2.11  reduces  to: 

Tipjb^o)^  [  (L,-  (V2+a)  L^)  a]  =2gbC^  (2.12) 

Substitute  K^.  =  bo)/U  for  w*  and  q  =  into 

equation  2.12  which  gives: 

2‘KgbKl(L^-(V2*a)L^)  a=2gbC^  (2.13) 

After  cancelling  and  siibstituting  Kp  for  K^.: 

[L„-(yz+a)Lja=C^  (2.14) 


This  relationship  can  be  further  broken  down  into  the  real  and 
imaginary  parts: 


I  mag: 


nK‘ 


[iL,-  (Yz+a)  iLj  a  =Q.„sin  (2.15) 


Plunge  Case:  of=0  using  equation  2.11  gives: 
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Real : 


(2.16) 


It 

^  [L,-  {V2^a)  LJ  a=C^.cos 

npb^<ji^Lf,—  =2gbCij^  (2.17) 

The  panel  code  uses  h/2b  for  analysis  because  it  uses  full 
chord  vice  half  chord.  Therefore  equation  2,17  becoines: 

=2(3i>C^  (2.18) 

Substituting  as  before  for  w: 

2ngbKlL^2  ( ^ )  =2gbC^  (2.19) 

Cancel  and  substitute  Kpi 

(  (  h  .  (2.20) 

This  can  also  be  broken  up  into  imaginary  and  real  parts  as 
before. 

MOMENT: 

Equating  equations  2.7  and  2.9  results  in: 

iipb«a)2(  [A4-(y2+a)Lj  |+[W.-(y2-*-a)  +  (y2+a) a)  ^2.21) 

=igb^yfR^ 

For  pitch:  h/b  =  0 
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Mp=4qb^yjR^^lle 


(2.22) 


(2.23) 


resulting  in: 

iipjb^o)2(  [Ar,-(y2+a)  (L„+M^) +(y2+a)2L^]  a)  =4gi)2C^  (2.24) 


After  substituting  and  cancelling: 


y>2 

^  [M^- (V2+a)  (L.+>f^)  +  (y2+a)  =C^  (2.25) 


REAL  Mh  =  ^ 


^  [M^- (V2*a)  (L.+y2)  +  {y2+a) 2lJ  =q«,cos (2*26) 
8 


IMAG:  Mh  =  0 


[iAf,-(y2+a)  (iLJ*{V2->-a)ULt,]=C^sin(^^)  (2.27) 

8 


For  plunge,  a  =  0,  equation  21  reduces  to  ; 


REAL;  Mjj  =  M 


-|K-p  ( ^)  [Ml,-  (y2+a)  L„]  =C^cos  (<l>«b) 


(2.28) 
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IHAG:  Mh  =  0 


(V2-^a)L^=C«,sin(V)  (2.29) 

Comparisons  are  shown  for  various  cases  of  pitch  and  plunge. 
The  tables  include  pitch  values  of  1  (Tables  2. 2 -2. 3)  and  6.7 
degrees  (Tables  2. 4-2. 5),  plunge  (h/2b)  values  of  .01  (Table 
2.8)  and  .0833  (Table  2. 6-2. 7).  The  graphs  include  1  degree 
pitch  (Figures  2.7-12,  2.13-20),  6.7  degree  pitch  (Figures 
2.21-28,2.29-36),  . 01  h/2b  plunge  (Figures  2 . 53 - 56 )  ,  and  .0833 
h/2b  plunge  (Figures  2.37-44,  2.45-52). 

4.  Results 

The  tables  and  graphs  show  that  the  panel  code 
predicts  the  Theodorsen  results  accurately.  An  initial 
question  that  was  first  addressed  for  the  comparisons  was  how 
many  cycles  to  use  for  good  consistent  phase  results.  Initial 
runs  were  made  at  several  different  cycle  values  and  the 
results  are  shown  in  Table  2.1.  It  can  be  seen  that  the  panel 
code  predicts  Theodorsen' s  results,  using  a  cycle  number  of 
three.  Increasing  the  cycle  number  just  takes  more  computer 
time  and  only  results  in  marginal  increases  in  accuracy. 

The  most  glaring  difference  appears  in  the  Ijj  Cjjjj 
comparisons  of  Figures  2.47-48.  The  panel  code  drops  off 
sharply  at  the  higher  end  of  the  reduced  frequency  spectrxim. 
The  reason  for  this  is  believed  to  be  due  to  the  magnitude  of 
h/2b  chosen  for  the  comparison.  The  code  was  rerun  for  a 
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Comparison  of  Phase  Calculations  Using  Various  Cycles. 

(Pitch,  6.7  deg.,  NACA  0007,  .37c,  SO  panels  top  and  bottom) 

Kp 

ItCyelis 

Cl  Phaso  Angle 

Cm  Phase  Angle 

Cl  Amplitude 

Cm  amplitui 

le 

1.00 

1 

182.0537 

-54.409 

1.00 

2 

208.1S92 

-46.424 

1.00 

3 

206.001 

•44.313 

0.51527 

0.083169 

%0Hf.2f3 

1.05% 

4.76% 

0.84% 

0.92% 

1.00 

4 

204.9365 

-43.33887 

0.51668 

0.082951 

1.00 

S 

204.3955 

-42.86817 

0.5172907 

0.082865 

_ 

1.00 

6 

204.0596 

-42.58497 

0.517598 

0.0828221 

1.00 

7 

203.8291 

-42.3916 

0.517776 

0.082792 

1.00 

8 

203.6031 

-42.25684 

0.51789 

0.082784 

%  OHf.  718 

Q.11% 

0.32% 

1 

3.60 

2 

264.44052 

59.79002 

1.0931 

wsm 

3 

261.9737 

58.29977 

1.09954 

0.21936 

%0iff.2/3 

0.94% 

2.56% 

0.59% 

3.35% 

mm 

260.17877 

57.4404 

1.10271 

0.21903 

1  3.60 

5 

259.2686 

56.79395 

1.104505 

0.218845 

wmm 

6 

258.5186 

56.3252 

1.105599 

0.218739 

3.60 

7 

257.9483 

•56.0098 

3.60 

8 

257.6163 

-55.380567 

1.106764 

0.218632 

%  0i(f.  7/8 

0.13% 

1.14% 

0.04% 

0.02% 

TABLE  2.1  PHASE  CALCULATION  VS  CYCLE  NUMBER 


series  of  h/2b  values  and  the  percent  difference  for  the  panel 
code  to  Theodorsen  was  plotted  in  Figure  2.4.  This  chart  shows 
that  the  h/2b  value  chosen  has  a  tremendous  impact  on  the  code 
results.  An  h/2b  value  of  .01  gave  an  acceptable  error  of  10% 
at  Kp  =  8.  Runs  were  completed  with  a  value  of  .01  h/2b  and 
the  favorable  results  are  shown  in  Figures  2.53-56. 
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o  o  o  o  o  o  o 

CsJ  O  00  <0  ^  CM 


U9SJOpOOl|X 
UIOJJ  leusd  k5  wa  % 


Figure  2.4  Effect  of  h/2b  values  on  In 
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RmI  Part  of  Cl  for  Panal  and  Thoodorton  (pitch,  1.0  dog.  .37c, 
NACA0007, 75  panels  top  and  bottom,  Scycloo  of  85  calc,  por  cyclo) 


TTiaodofsta 


/  ,1 


i  ■ 

! 

2  - - 

1 

i 

1 

1 

0 

0.5 

1.5  2  2.5 

3 

3.5  i 

Xp  -  2MU 

Figure  2.6  1  degree  pitch  C,  Re 
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29 


Imag.  Part  of  Ci  for  Paaal  aad  Thaordoraan  (pitch,  1.0  dag,  .37e, 
NACA0007, 7S  panala  tap  and  bottom,  Seydos  of  85  calc,  par  eycia) 
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Mag.  of  Cl  for  Paaal  and  Thaordorsan  (pitch.  1.0  dag,  .37c, 
NACA0007. 75  panaia  top  and  bottom,  Scyclac  of  65  calc,  par  cycio) 


CL  Moq. 


Mag.  of  Cl  for  Panel  and  Thaordorsan  (pitch,  1.0  dag,  .37e, 
NACA0007, 75  panala  top  and  bottom.  Scyeloc  of  65  calc,  par  cycle) 


Figure  2.10  1  Degree  pitch  Cj,  Magnitude 
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1  Degree  pitch  Phase 


iml  and  ThaerdorMn  (pitch,  1.0  dag,  .37c, 

I  tap  and  bottam,  Scyclas  of  85  calc,  par  cycia) 


p 

fnelCode 

_ z 

1 

■ 

/ 

r 

Theodor. 

en 

1  1.5  2  2.5 

Ka  -  2bw|U 


Figure  2.12  1  degree  pitch  Cj,  Phase 
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35 


(pitch.  1.0  deg.  .37c. 
of  65  calc,  per  cycle) 


1 


Cm  hnag. 


Cm  Imag  vs  Kp  for  Panel  and  Theodorsen  (pitch,  1 .0 
deg.  .37c,  NACA0007.  75  panels  top  and  bottom, 
3cycles  of  65  calc,  per  cycle) 

Kp  =  2bw/U 


0  0.5  1  1.5  2  2.5 


Figure  2.16  1  Degree  pitch  Ck  Im 
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Cm  Phase 


Cm  Phase 


Figure  2.20  1  Degree  pitch  C,,  Phase 
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Cl  Reol 


Ci  Real  vs  Kp  Panel  Cede  (Pitch,  6.7  dag,  .37c,  NACA  0007, 50  panels 

tap  and  battam) 


Figure  2.21  6.7  degrees  pitch  Re 


Cl  Rul  vs  Kp  Pans!  Cods  (Pitch.  6.7  dog,  .37,  NACA  0007. 75  panels 

top  and  bottom) 


Figure  2.22  6.7  degrees  pitch  Re 


Cl  Imog. 


Cl  Imag.  vs  Kp  Paaal  Coda  (pheb.  8.7  dag,  .37g,  NACA  0007, 75  pansi 

tap  aad  batton) 


0  0.2  0.4  0.6  0.8  1  1.2  1.4  1.6  1.8  2  2.2  2.4 

Kp  •  2bw/U 


Figure  2.24  6.7  degrees  pitch  Im 
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Cl  Magnitude 


Figure  2.25  6.7  degrees  pitch  Cj,  Magnitude 
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53 


Cm  Imaginary 


57 


Cm  Magnhuda 


Cm  Magnitude  vs  Kp  Panel  Code  (pitch,  6.7deg,  .37c, 
NACA  0007,  50  panels  top  and  bottom) 


Figure  2.33  6.7  Degrees  pitch  Cj,  magnitude 
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Cm  MagnHucto 


Figure  2.34  6.7  Degrees  pitch  magnitude 
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Cm  Phas* 


Cm  Phase  vs  Kp  Panel  Code  (pitch,  6.7deg,  .37c, 
NACA  0007,  SO  panels  top  and  bottom) 


Cm  Phase 


Cm  Phase  vs  Kp  Panel  Code  (pitch,  6.7deg,  .37c, 
NACA  0007,  50  panels  top  and  bottom) 


Figure  2.36  6.7  Degrees  pitch  phase 
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62 


CL  RmI 


Real  Part  of  Cl  for  Panel  and  Theodorsen  (plunge, 
.0833  h/2b,  .37c,  NACA0007,75  panels  top  and 
bottom,  3cycle8  of  65  calculations) 


Figure  2.37  Pliinge  h/2b>.0833  Re 


CL  RmI 


Real  Part  of  Cl  for  Panel  and  Theodoraen  (plunge, 
.0833  h/2b,  .37c,  NACA0007,75  panels  top  and 
bottom,  3cycle8  of  65  calculations) 


Figure  2.38  Plunge  h/2b-.0833  Re 
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Imag  Part  of  Cl  for  Panel  and  Thaodorsan  (plunge, 
.0833  h/2b.  .37c,  NACA0007,75  panels  top  and 
bottom,  3cyclea  of  65  calctdations) 


Figure  2.39  Plunge  h/2ba.0833  C,  Im 

If 
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Imag  Part  of  Cl  for  Panel  and  Theodorsen  (plunge, 
.0833  h/2b,  .37c,  NACA0007,75  panels  top  and 
tX)ttom,  3cycles  of  65  calculations) 


0 


-0.1 


-0.2 


-0.5 


-0.6 


-0.7 


0  0.5  1  1.5  2  2.5 


Kp  =  2bw/U 


Figure  2.40  Plunge  h/2b-.0833  C,  Im 
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Magnitude  of  d  for  Panel  and  Theodorsen  (plunge. 
.0833  h/2b,  .37c.  NACA0007.75  panels  top  and 
bottom,  3cycle8  of  65  calculations) 
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Magnitude  of  D  for  Panel  and  Theodorsen  (plunge, 
.0833  h/2b.  .37c,  NACA0007,75  panels  top  and 
bottom,  3cycles  of  65  calculations) 


Kp  =  2bw/U 


Figure  2.42  Plimge  h/2ba.0833  magnitude 


CL  Phase 
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CL  Phase 


Figure  2.45  Plunge  h/2b-.0833  Re 


CM  Real 


Real  Part  of  CM  for 
.0833  h/2b,  .37c, 
bottom. 


! 

! 

I 

j 

1 - 

Figure  2.46  Plunge  h/2b>.0833  Re 


(•]  11*1 


73 


Imag.  Part  of  CM  for  Panel  and  Theodorsen  (plunge. 
.0833  h/2b.  .37c,  NACA0007,75  panels  top  and 
bottom,  3cycle8  of  65  calculations) 


CM  Imag. 


Magnitude  of  CM  for  Panel  and  Theodorsen  (plunge, 
.0833  h/2b,  .37c,  NACA0007,75  panels  top  and 
bottom,  3cyde8  of  65  calculations) 


Figure  2.49  Plunge  h/2ba.0833  C,,  Magnitude 


CM  MAG. 


Magnitude  of  CM  for  Panel  and  Theodoraen  (plunge. 
.0833  li/2b,  .37c,  NACA0007,75  panels  top  and 
bottom.  3cycles  of  65  calculations) 


0  0.5  1  1.5  2 

Kp  =  2bw/U 


i 


I 
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Figure  2.50  Plunge  h/2b>.0833  Magnitude 


Phase  of  CM  for  Panel  and  Theodorsen  (plunge, 
.0833  h/2b,  .37c,  NACA0007,75  panels  top  and 
bottom,  3cycles  of  65  calculations) 


Figure  2.51  Plunge  h/2b-.0833  C,,  Phase 
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CM  Phase 


Phase  of  CM  for  Panel  and  Theodorsen  (plunge, 
.0833  h/2b,  .37c.  NACA(X)07,75  panels  top  and 
bottom,  3cyclas  of  65  calculations) 


Kp  2bw/U 


Figure  2.52  Plunge  h/2ba.0833  C,{  Phase 
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80 


82 


.55  Plimge  h/2bB.01  Magnitude 
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Phase  of  CM  for  Panel  and  Theodoraen  (plunge,  .01 
h/2b,  .37c,  NACA0007,100  panels  top  and  bottom, 
Scydes  of  65  calculations) 


Figure  2.56  Plunge  h/2b>.01  Cj,  Phase 


Comparison  of  CM  Imag  vs  Kp  for  Various  Valuos  of  Pitch 
(NACA0007.100  pansls,1^and  6.7  dag.  pitch.3cye65M0c,) 
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In.  My 


'Comparison  of  Imaginary  Dimansionlots  Aarodynamic  Coefficient  dm 
My»8*lm  Cmalpha|(alpha*pi*Kp*2)(NACA0007,10Q  panals.l. 3,6.7  dag. 

piteli,3cyc65,0c.) 


80  - 


Kp  -  2bw/U 


Figure  2.58  Dimensionless  Aerodynamic  Coefficient  for 
1 . 0 , 3 . 0 , and  6 . 7  degrees 
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Ill .  FLIITTER  DETERMINANT 


The  proven  accuracy  of  the  UPOT  Code  enabled  it  to  be  used 
for  the  solution  of  the  flutter  determinant. 

A.  FLUTTER  THEORY 

In  order  to  analyze  the  phenomenon  of  flutter,  it  is 
necessary  to  obtain  the  equations  of  motion  of  the  system.  To 
simplify  the  problem  the  assumption  is  made  that  the  actual 
motion  of  the  system  can  be  considered  a  combination  of 
fundamental  wing  bending,  and  fundamental  wing  torsion.  The 
system  can  then  be  replaced  by  an  equivalent  system  containing 
an  airfoil  section  of  unit  span  restrained  by  springs  against 
independent  vertical  motion  (bending) ,and  torsion  as 
illustrated  in  Figure  3.1.  This  paper  will  not  consider  the 
aileron  hinge  case  so  S  and  c  are  set  equal  to  zero.  According 
to  the  class  notes  of  M.  Platzer  [ref.l]  the  formulation 
proceeds  as  follows: 

Consider  the  balance  of  the  elastic,  inertial  and 
aerodynamic  forces  on  a  mass  element: 

•  Total  Inertial  force:  -/ dm(h"+rQ'" )  =- (Mh"+S„a" ) 

Mass :  M=/  dm 

Static  Moment  about  the  elastic  axis:  =  I  rdm 

•  The  moments  about  the  elastic  axis  are: 
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-Ir(h"+ra")dw  =  - 


(3.1) 


Mass  moment  of  inertia  about  elastic  aocis 
I„  =  Jr^dm 


•  Elastic  restoring  forces  are:  -hCjj  "QfCa 
The  Equations  of  motion  therefore  become: 


=  L 

a"l^-^h"s^+aC^  =  M 


(3.2) 


Where:  =  Torsional  stiffness  of  the  wing 

Cjj  =  Stiffness  of  the  wing  in  translation  (plunge) 
M  =  Mass  of  the  wing  per  unit  span 
These  equations  can  be  written  in  a  different  way  by 
expressing  the  spring  constants  in  terms  of  the  natural 
frequencies.  Consider  the  airfoil  to  be  so  restrained  that 
only  one  degree  of  freedom  is  permitted.  The  equations  of 
motion  become: 


Mh"+hC^  =  0 


so  that 


=  0 


so  that 


(3.3) 

1 


Hence:  Cj,  =  Moj^^ 

The  small  structural  daunping  of  metal  aircraft  may  be 
approximated  by  a  force  that  opposes  the  motion  and  is  in 
phase  with  the  velocitv.  One  assumes  therefore  that  the 
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magnitude  of  the  damping  is  proportional  to  the  elastic 
restoring  force.  Since  the  motion  of  the  airfoil  is  harmonic 
at  the  critical  flutter  condition,  the  structural  damping  can 
be  accounted  for  by  replacing  the  terms: 

hC^  with  hC^d+igj,) 
with  aC^d+ig^j) 

Where  g^  and  g^  are  damping  constants  multiplied  by  i  to  ensure 
that  the  damping  force  is  in  phase  with  the  velocities  in  the 
simple  harmonic  motion. 

From  equation  3.3  with: 

h(t)  =  he'“*  and  a(t)  =  ore’"^ 

We  have: 

h"  =  <i)^he'“*  and  a"  = 

And  the  equations  of  motion  become: 

=L  (3.4) 

( -w^a  J,-w^/35,+oC„)  =  M  (3.5) 

The  equations  for  the  aerodynamic  forces  were  given  by 
Fung  [ref. 5]  and  are  shown  here; 

L=npb^<^^  +  [L.-  ( 1  +a)  LJ  a)  (3.6) 

Equating  equation  3,4  to  3.6  and  3.5  to  3.7  yields: 

Substituting  into  equation  3.8  and  3.9  for  C,,  and  and 
using  the  following  dimensional  terms; 


89 


(3.7) 


M=npb*<o^  (  [M^,-{  —  +a)  Lj  ^  + 

[Af„-  (-i+a)  iL^+M„)  +  (i+a)  a)e^^‘= 

(-<ji^Mfi-<>}^aS„+hCjj)  =i:pJb^ci^  (L}^-^+  [L^- (--+3)  a) 

JD  A 

(-w2aI,-o>2ii5,+oC„)  =npi)^w2(  [M^- ( -l+a)  LJ  ^  + 
[M.-(-|+a)  (L„+M*)+(^+a)2Lja) 


(3.8) 


(3.9) 


H  = 


Af 

tipjb^ 


(3.10) 


The  equations  simplify  and  after  bringing  all  terms  to  the 
left  result  in: 

A^+Ba=0 

b 

(3.11) 

D-?  +£a=0 


This  is  a  homogeneous  equation  whose  solution  is  obtained 
if  the  flutter  determinant  is  zero. 


A 

D 


(3.12) 


Where : 
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(3.13) 


(O  0),  "  " 

S=n;f„+L.-L^(Va+a) 

D=\ix„+V2-Lj,  iVz+a) 

E=\ir^^  [1-  ( — )  (l+ig„)  ]  -V2  (Vs+a)  ('h+a)  *Lf^  (Va+a)  ^ 

0) 

/i  is  the  ratio  of  the  mass  of  the  wing  to  the  mass  of  a 
cylinder  of  air  of  a  diameter  equal  to  the  chord  of  the  wing. 

and  Ujj  are  the  natural  angular  frequency  (rad/sec)  of 
torsional  vibration  around  "a"  (elastic  axis)  and  the  natural 
frequency  in  deflection,  respectively,  is  the  location  of 
center  of  gravity  of  the  wing  measured  from  a.  w  is  the 
circular  frequency  of  wing  vibration. 

The  relationships  between  the  code  and  Theodorsen  derived 
earlier  in  Chapter  II  can  be  used  here  to  simplify  the 
equations:  (note:  no  damping  in  this  case  g„=gf,=0) 

For  A:  manipulating  equation  2.20, 


TTh 


(3.14) 


resulting  in; 


(1)  (<>_ 


2C, 


Ls 


(3.15) 


For  B:  manipulating  equation  2.14, 
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L.-L*(y2+a) 


_ 

TlKpU 


resulting  in: 


B=lix,+  ^ 


fli  o 

IlKpO. 


For  D:  manipulating  equation  2.28 


i„-L^(y2+a)  =. 


4C. 


Mh 


(3.16) 


(3.17) 


(3.18) 


resulting  in: 


D=\ix, 


4C. 


Mh 


(3.19) 


For  E:  manipulating  equation  2.24 


8  r? 

-1/2  (1/2+3)  (V2+a)  {Vz+a)  - — 

anKp^ 


(3.20) 


resulting  in: 

E=nri  [1- ( -^ ) ^  (3.21) 

w  anKp^ 

The  determinant  is  expanded  to  AE-BD=0,  and  the  real  and 
imaginary  parts  are  set  equal  to  zero.  Substituting  (03^/03)^=^ 
X  and  solving  the  real  (2  roots)  and  imaginary  (1  root) 
equations  for  values  of  X  corresponding  to  each  reduced 
frequency  value.  These  X  values  can  be  plotted  as  SQRT{X) 
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against  Kp  and  any  intersections  of  real  and  imaginary  parts 
signify  a  flutter  point. 

Knowing  that  : 


,  _  2jb<o 
U 


and  <JX=  ( ) 

(0 


solve  for  U 


critical 


U. 


cri  tical 


which  is  the  critical  flutter  speed. 


(3.22) 


(3.23) 


B .  UPOTPLDT  CODE 

1.  FORMULATION  AND  INPUT 

The  equations  derived  in  the  flutter  theory  section 
above  were  programmed  into  a  FORTRAN  subroutine  and  attached 
to  the  UPOT.f  code.  The  UPOT  code  was  modified  first  to 
enable  it  to  conduct  a  frequency  sweep  of  pitch  and  plunge 
simultaneously.  The  resulting  frequency  sweep  pitch  and 
plunge  array  data  is  then  sent  to  the  flutter  subroutine  which 
provides  the  values  of  SQRT{x)  and  Kp  for  plotting.  The 
program  also  gives  a  best  guess  for  the  based  on  the 

difference  betweeen  the  real  and  imaginary  SQRT  (X)  values. 
The  input  file  UPOTFLUT.IN  is  very  similar  to  the  regular 
UPOT. IN  file  with  the  addition  of  actual  physical  properties 
of  the  system  being  analyzed.  The  user  should  start  the 
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analysis  in  the  pitch  mode  first  (lOSCIL  =1,  ITRANS  =0)  to 
ensure  complete  coverage  of  all  frequencies  of  interest.  The 
following  relations  were  taken  from  NACA  TR-685  [ref. 8]  and 
should  prove  helpful  in  determining  the  physical  properties 
needed  for  program  operation. 

K  =  mass  ratio  =  irpb^/M 

K  =  1//I 

x„  =  S„/Mb 

r„2  =  i^/Mb2 

2 .  OUTPUT 

Outputs  from  the  code  have  been  limited  to  reduce  the 
amount  of  computer  space  taken  up  by  the  code  operation.  A 
scimple  input  and  output  file  are  contained  in  Figure  3.2. 
The  following  list  describes  the  input/output  files  and  the 
data  they  contain: 

a.  UPOTFLUT.IN  The  input  file  figure  3.2a 

b.  CL.d  Scune  as  UPOT.f  output 

c.  CM.d  Same  as  UPOT.f  output 

d.  PHASE. d  Same  as  UPOT.f  output 

e.  CPSS.d  Same  as  UPOT.f  output 

f.  CPUOOB.d  Soune  as  UPOT.f  output 

g.  PHZSWP.d  This  file  contains  Kp, ^La* 

h.  PLHZSWP.d  This  file  contains  Kp,  ^Lh' ^Mh 

i.  PITCH. in  This  file  contains  Kp  Re,  Im,Cjjjj,  Re, 

CmoN  im 
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b  =  semicliord  (ft.) 

cb  —  distance  between  midcliord  and  aileron  hinge,  positive  if  aft  of 
midcliord  (ft.) 

eb  distance  between  midchord  and  aileron  leading  edge,  positive 
•  aft  of  midcliord  (ft.) 

^cib  >=  distance  between  rotation  point  (elastic  axis)  and  midchord, 
positive  if  aft  of  midchord  (ft.) 

h  —  bending  deflection  of  rotation  point  (elastic  axis),  positive  down¬ 
ward  (ft.) 

( a  angular  deflection  about  rotation  point  (elastic  axis),  positive 
for  leading  edge  up  (radians) 

0  <=  angular  deflection  of  aileron  about  aileron  hinge  relative  to  wing 
chord,  positive  for  aileron  leading  edge  up  (r<idiuns) 


Figure  3.1  Slng>ll£led  System  Geometry 
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5^  4 


itdin 


AIRFOIL  TYPE  ;  NACA  0012  AIRFOIL 
NLOWER  V  50  ,  NUPPCR  >  SO 


IFl.AG 

NLOWER 

NUPPER 

0 

50 

50 

A^PFOIL  TYPE 

i.^P 

lOSCTL 

ALPl 

ALPMAX 

PIVOT 

0 

1 

-3,0 

3.0 

0.3 

PREO 

RFOSTP 

RFQFNL 

.85 

.01 

.95 

IGUST 

UGUST 

VCUST 

0 

0. 

0. 

TTRANS 

OELMX 

DCLHY 

DELI  PHASE 

0 

0.00 

.0833  - 

.0833  0.00 

CYCLE 

NTCYCLF 

TOI. 

3 

85 

0,005 

naot  t 

naot  X 

aoa  values  multiplied 

by  10  (Integer) 

2  05  10  20  25  39  50 

S«m{ -chord  Walpha  Hh  Mass 

«  90  72.0  .53799 

lalpha  Salpha  Density 

4.84102  .6A546B  .002378 

Convnent  s  .  .  . 


TRAMP  Ot  n/a 

1:  Straight  ramp 
2:  Modified  ramp 


RFRCO  Is  based  on  fu))  chord 


TOSCIL  0:  n/a  RfRCQ  Is  based  on  iuH  chord 

1:  Sinusoidal  pitch,  motion  starts  at  min  Ana 

7TRANS  0:  n/a 

1:  Trans lat  lonal  harmonic  oscillation 

CYCLE  :  f  of  cycles  for  oscillatory  motions 

-In  case  of  ramp,  cycle-l«S  denotas  airfoil  is  held 
at  max  aoa  for  the  duration  of  .5  cycle 
-For  steady  state  solution  set  It  to  0 

NTCYCLE:  f  of  time  steps  for  each  cycle 

CYCLE*NTCYCLE  Is  limited  to  200  currently. 

NAOT:  I  of  Input  aoa  for  cp  output 

-  angles  should  be  In  Increasing  order, 

-  for  oscllatory  motions  angles  should  increase 
first,  then  decrease.  Decreasing  angles  are  for 
the  return  cycle., 

SEMI-tfHOFD  H*lf  Chord  In  foot , 

Mslpha,Wh,  uncoupled  natural  frequencies  of  the  system  in  question. 
Malpha  is  pitch  and  Hh  Is  plunge(HZ). 

Mass  specific  mass  of  the  system  In  slugs/foot  of  span 

lalpha  Moment  of  Inertia  of  system  about  the  elastic  axis  (a) 

per  unit  span  length. 

Salpha  Static  moment  of  wlng-alleron  per  unit  span  length 

Density  Mass  of  air  per  unit  of  volumelslugs  per  ft*3) 


Figure  3.2a  UPOTFLUT . in  exaiig>le  Input  file 
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Figure  3.2b  UPOTPLUT  exas^le  output  file 
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. . .  "  " 


'  uf,,  ^  Papa  jj 


kp<  0.S700000  Ifreq  3 

AMPLITUDE;  clamp,  cmamp  :  0.J11T9JS  2.S7655HE-0? 

loscll  >  Oltrana  -  0 


PHASE;  clp,  cnp 

AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  HBAR 


203.3604  '64.02047 

1.5ai342BE-03  0.1048S66 


-0.1365460 


-1.163S614E-02 


PHASE  SHIFT  ANALYSIS 
FREO  *  0.8700000 

w  0.8700000 

kp°  0.8700000  Ifreq  3 


AMPLITUDE;  clamp,  cmamp  : 
loscll  «  Oitrans 
PHASE:  clp,  cnp  : 
AVERAGE  DRAG,  TOTAL  DRAG  ; 
ETAS,  NBAR  : 


0.2830637  2.2B74046E-02 


1 

271.0655  36.47753 

-6.4334869E-03  -0.4246101 
0.1045149  -0.1231114 


FREO  SWEEP 
FREO  *  0.880000 


PHASE  SHIFT  ANALYSIS 
FREO  >  0.8800000 

w  0.8800000 

kp«  0.8800000  Ifreq  4 

AMPLITUDE;  clamp,  cmamp  :  0.2318744  3.0O371B0E-02 

loscll  '  Oitrans  =  0 

PHASE;  clp,  cmp  :  203.8330  -64.15525 

AVERAGE  DRAG,  TOTAL  DRAG  ;  1.58804258-03  0.1048108 

ETAS,  NBAR  :  -0.1335113  -1 . 1894441E-02 


PHASE  SHIFT  ANALYSIS 
FREO  •  0.8800000 

w  0.8800000 

kp>  0.8800000  Ifreq  4 


AMPLITUDE;  clamp,  cmamp  ; 
loscll  •  Oitrans 
PHASE;  clp,  cmp  : 
AVERAGE  DRAG,  TOTAL  DRAG  ; 
ETAS,  NBAR  : 


0.2856959 
•  1 
271.3546 
-6.54S8103E-03 
0.1042214 


2.3258235E-02 

36.17284 

-0.4322215 

-0.1256711 


FREO  SWEEP 
FREO  ■  0.890000 


PHASE  SHIFT  ANALYSIS 


FREO  - 

H  0.8900000 

kp<  0.8900000  Ifreq 

AMPLITUDE;  clamp,  cmamp  : 
loscll  ■■  Oitrans 

PHASE;  clp,  cmp  : 

AVERAGE  DRAG,  TOTAL  DRAG  : 
ETAS,  NBAR  : 


i.  8900000 

5 

0.2319494 
-  0 

3.0311935E-02 

204.3115 

-64.29196 

1.5859993E-03 

0.1046759 

-0.1304959 

-1.2153637E-02 

PHASE  SHIFT  ANALYSIS 
'  FREO  •  0.8900000 

w  0.8900000 

kp-  0.8900000  Ifreq  5 

AMPLITUDE;  clan«>,  cmamp  ;  0.2881391  2 . 3644408E-02 


Flgur«  3.?c  UPOIFLUT  •xainpl«  output  £11« 
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^  ;  I 


loscll  >  Oltrans 

PHASE;  clp,  cmp 

AVERAGE  DRAG,  TOTAL  DRAG 
ETAS,  MBAR 


FREQ  SWEEP 
FREQ  »  0.900000 


271.6436  35.87401 

'6.6649993E-03  -0.4398900 


0.1039310 


-0.1282581 


w  0.9000000 
kp=  0.9000000 


PHASE  SHIFT  ANALYSIS 
FREQ  ^  0,9000000 

Ifreq  6 


AMPLITUDE;  clamp,  cmamp  :  0.2320280 

lo9cil  =  Oltrans  =■  0 

PHASE;  dp,  cmp  204.7803 

AVERAGE  DRAG,  TOTAL  DRAG  :  1 .Se2e811E-03  0.1044701 


3.0592278e-02 
-64.43259 


ETAS,  W8AR 


-0. 1274617 


-1.2418482P.-02 


PHASE  SHIFT  ANALYSIS 
FREQ  r  0.9000000 


I  £r«q 


w  0.9000000 
fcp-  0.9000000 

AMPLITUDE;  damp,  cmamp  :  0.2909827  2.4039967E-02 

load!  »  Oltrans  s  1 

PHASE;  dp,  cmp  :  271.9288  35.58691 

AVERAGE  DRAG,  TOTAL  DRAG  :  -6 . 782107ie-03  -0.4476191 

ETAS,  WBAR  :  0.1036481  -0.1308680 


FREQ  SWEEP 
FREQ  s  0.910000 


PHASE  SHIFT 

ANALYSIS 

FREQ  ^  0 

,9100000 

H  0.9100000 
kp.  0.9100000 

1  f  req 

7 

AMPLITUDE;  damp,  cmamp  : 

0.2321466 

3.0865142E-02 

loacll  « 

Olt  rans 

*  0 

PHASE;  dp, 

cmp  : 

205.2529 

-64.61227 

AVERAGE  DRAG,  TOTAL  DRAG  : 

1 .5800473E-03 

0.1042831 

ETAS,  WBAR 

* 

-0.1245367 

-1.2687406E-07 

* 

PHASE  SHIFT 

ANALYSIS 

FREQ  *  0 

,9100000 

w  0.9100000 
kp=  0,9100000 

1  f  req 

7 

AMPLITUDE;  damp,  cmamp  : 

0.2936268 

2.4440434E-02 

loacll  ■ 

Olttans 

=  1 

PHASE;  dp. 

cmp  : 

272.2159 

35.29980 

AVERAGE  DRAG,  TOTAL  DRAG  : 

-6.900007ee-Q3 

-0.4554005 

ETAS,  WBAR 

: 

0,1033707 

-0.1335003 

FREQ  SWEEP 
FREQ  »  0.920000 


PHASE  SHIFT  ANALYSIS 
*  FREQ  «  0.9200000 

w  0.9200000 

kp*  0.9200000  l£feq  8 


Figure  3. 2d  UPOTFLUT  exan^le  output  £lle 
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W»rtir>Or.:  Cl-imp,  cmii^  :  S.11?ift^7r-0? 

loseil  *  OltruAs  ~  0 

rRASCi  elp,  cmp  205.1313  -ft4.(9(2i 

AVtRAOe  DRAO.  TOTM.  ORAO  t  1 .51i2 3  >2r“0.t  0.1041<37 

r.TA5.  WRAP  t  -0.1HRU3  .1.2MM*ir-02 


rpAse  sprrr  awalmis 

fpr.O  *  0.9200000 

w  0.1200000 

kp-  0.9200000  • 

AMFLITOOr.:  e««»p  :  0.291^2105  2 . 4P415§1F-02 

io«cll  ”  Oitrjtns  *  1 

rHASE;  dp,  cmp  2i2.50tl  35, 01200 

AVCRAOE  OAAC,  TOTAI  ORAO  :  '1.0100030F>01  •0.4(32331 

riAS,  WHAP  0.1030901  -O.M(154( 


rpEO  swrrp 

rpro  -  0,910000 


rHASE  SKirt  AKALVSIS 
IPfO  -  0.9299999 

w  0.9299999 

kp^  0.9299999  9 

AWLiTtiOr..*  clump,  e«t«mp  :  0.212959*  1.M11310r-02 

ioccit  *  01tr«n»  -  D 

fRASCl  etp,  cap  20(.2090  -(9.12113 

AVEpAoe  OMC,  TOTAL  0AA<?  s  1 . 5155(24r-03  0.1039011 

r.TA5.  W*A*  t  -0.1191(51  .|.3221(01E-02 


PHASE  Sltirt  AWAUSIS 
fprO  •  0.9299999 

w  0.9299999 

fcp-r  0.9299999  9 

AHTLItOOn  diHi»p,  CiMPp  t  0.2909111  2.52901122-02 

loaell  -  OUrPAt  *  \ 

PHASE;.  cfp,  rap  ;  212.2901 

AVEAAOC  OHAO,  TOTAL  OPAO  j  -1.1302109r.-03  -0.4)11219 

ETAS.  WPA*  t  0.)02*32«  -O.nSilM 


rOEO  SWEfP 
rpr.Q  <  0,940000 


PHASE  SNtrt  AIIALYSI9 
fpro  -  0.9199999 

a  0.9399999 

kp-  0.9199999  lfrr<!  10 

MirLlTUOf;  d«mp,  capap  j  0.2324414  1.  )(990H1:-02 

iopcll  »■  OUtPAP  •  O 

PMSC;  dp.  20(.(94(  -(4,94021 

AVERACC  drag,  total  PHAC  t  1.51l5122r-03  0.1031211 

ETAS.  WHA*  J  -0.11(415*  -1.14949P9E-02 


PRASE  SRIfT  AWALTStS 
rRf.O  0.9599999 

a  0.9199999 

kp^  0.9199999  Ifrrn  10 

AMPtnoor.}  rlA«p,  empap  :  0.1015(90 

iQfldl  «  OttCPAS  «  I 


2  5(599(?r-02 


PHASE;  rip,  rap  t  213.0)9 

AVERAoe  ORAr,  TOIAI,  npAC  :  -1.250)1 

Tins,  WPA*  0.1025) 

Httabac  of  Kp  Valur«  • 
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Figure  3.2e  UPOTFLOT  example  output  file 
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j.  PLUNGE. in  This  file  contains  K_,  C^h  RE,  C^h  IM,  RE, 

CMh  IM 

k.  FLUTPLOT.d  This  file  contains  IC,SQRT(x)  Re,  SQRT(X)  RE, 

SQRT(X)  IM 

3 .  VALIDATION 

The  program  was  tested  against  some  sample  cases  to 
check  for  code  validity.  The  first  case  was  taken  from 
reference  6  excimple  #1,  p.  236.  Figures  3.3  and  3.4  show 
plots  of  the  FLUTPOT.d  file.  Figure  3.3  shows  the  initial 
look  over  a  wide  range  of  Kp  and  after  finding  the  approximate 
flutter  location  Figure  3.4  shows  a  closer  look  at  the  Kp 
range  of  interest.  This  example  calculated  a  of 

161.985  ft/sec.  which  compares  favorably  to  the  example  value 
of  162  ft/sec.  The  next  example  was  taken  from  NACA  TR-685 
[ref. 8]  case  #1  p.  8.  Figures  3.5  and  3.6  again  show  the 
initial  and  final  looks  for  this  analysis.  The  example  called 
for  a  of  567  miles/hr  and  the  program  returned  a 

value  of  570  miles/hr.  Next,  the  code  was  tested  over  a  range 
of  i^atios  as  done  in  NACA  TR-685,  pll.,  graph  I-A(a)  . 

Figure  3 . 7  shows  the  comparison  between  the  two  methods . 
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Figure  3.4  Final  look  at  flutter,  exas^le  1 
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Kp=2bwAJ 


Kpanel  =  2bwAJ 


Flutter  Ex.  (N0007,' 


Figure  3.6  Final  look  at  flutter,  example  2 
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Kpanel  =  2bw/U 


TV.  FLOW  VISUALIZATION  EXPERIMENT 


A.  INTRODUCTION 

The  purpose  of  this  experiment  was  to  document  the 
production  of  thrust  by  a  plunging  airfoil.  This  was  a 
preliminary  experiment  to  better  understand  the  vortex  pattern 
produced  by  a  plunging  airfoil,  and  to  examine  the  production 
of  thrust  using  smoke  flow  visualization  techniques. 

An  explanation  of  what  constitutes  a  propulsive  vortical 
signature  along  with  smoke  flow  visualization  of  the 
propulsive  vortical  patterns  is  given  in  Reference  7.  In  this 
reference,  the  explanation  is  given  by  contrasting  the 
vortical  pattern  produced  by  a  cylinder  (drag)  with  the 
vortical  pattern  produced  by  a  plunging  airfoil  (thrust)  .  The 
cylinder  produced  a  vortical  sheet  where  the  top  row  of 
vortices  rotated  clockwise  and  the  bottom  row  of  vortices 
rotated  counterclockwise.  This  pattern  induces  a  velocity 
component  in  the  upstream  direction  (Biot-Savart  law) .  In 
contrast,  the  plunging  airfoil  produced  a  clockwise  rotating 
vortex  sheet  on  the  bottom  row.  This  pattern  induces  a 
velocity  component  in  the  downstream  direction.  Reproduction 
of  the  flow  visualization  data  from  Reference  7  is  shown  in 
Figure  4.1. 
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Figure  4.1  View  of  flow  over  cylinder 
airfoil  (bottom)  [Ref.  7] 


B.  THEORY 


A  comparison  was  done  using  the  incompressible  panel  code, 
U2DIIF.  The  purpose  of  this  study  was  to  examine  the  vortical 
pattern  produced  by  the  panel  code,  and  determine  if  the 
vortical  signature  matched  experimental  results.  The  input  to 
the  panel  code  was  set  up  to  best  match  the  conditions  of  the 
experiment  described  in  the  next  section.  The  panel  code  was 
run  using  a  plunge  amplitude,  h/2b  equal  to  .1977,  a  reduced 
frequency  of  1.8  and  a  zero  mean  AOA.  The  results  of  the 
vortical  pattern  are  shown  in  Figure  4.2.  Aside  from  the 
starting  vortex,  this  is  clearly  a  thrust  producing  vortical 
sheet.  Furthermore,  the  vortical  pattern  is  similar  to  that 
produced  by  the  experiment  shown  in  Figures  4.10  and  4.11. 

C.  EXPERIMENTAL  SETUP 
1.  Plunging  Airfoil 

The  plunging  airfoil  used  in  this  paper  was  originally 
a  wing  taken  from  the  rotor  of  a  model  helicopter.  The  wing 
was  attached  to  a  MB250  Shaker  Table  as  shown  in  Figure  4.3 
The  wing  was  made  from  a  NACA0007  airfoil  section  and 
consisted  of  a  2.45"  chord  and  a  22"  span.  The  wing  was  built 
from  a  foam  core  and  finished  with  a  layer  of  graphite  epoxy 
composite  for  added  fatigue  strength.  The  airfoil's  drive 
mechanism  was  a  MB  250  Shaker  Table  capable  of  1"  total 
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deflection.  The  Shaker  Table  was  limited  only  by  resonance 
frequencies  of  the  wing  which  occurred  around  20  Hz  or  1200 
rpm. 

2 .  WIND  TDMMEL 

The  wind  tunnel  used  in  this  experiment  was  a  very  low 
speed,  low  turbulence  smoke  tunnel.  It  is  made  of  plexiglass 
walls  and  a  contraction  ratio  of  2.8:1.  The  motor  provides 
wind  tunnel  velocities  between  0  and  10  feet  per  second  (fps)  . 
The  smoke  was  created  using  a  Rosco  smoke  generator  and  piped 
into  the  tunnel  in  the  test  section  using  a  small  seven  tube 
smoke  rake  constructed  for  this  experiment.  Figure  4.4  is  a 
photograph  of  the  wind  tunnel  and  smoke  rake  used  in  this 
experiment . 

D.  TEST  PROCEDURE 

Testing  was  conducted  in  the  low  speed  smoke  tunnel  under 
several  different  conditions.  The  speeds  of  the  tunnel  were 
approximately  1.04  fps,  1.47  fps,  and  1.56  fps  (measured 
visually) .  These  low  speeds  allowed  good  pictures  and  the 
cdDility  to  get  higher  reduced  frequencies  without  calling  for 
too  high  a  load  on  the  wing.  The  actual  plunging  harmonic 
frequencies  ranged  from  l  to  15  Hz  and  amplitudes  from  1/16" 
to  1"  peak  to  peak.  The  tunnel  was  initially  turned  off  and 
the  Shaker  Tadale  turned  on  with  stagnant  smoke  in  the  tunnel . 
The  purpose  was  to  see  if  the 

plunging  airfoil  would  draw  the  smoke  through  the  tunnel  like 
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a  fan,  thus  showing  the  production  of  thrust  by  the  plunging 
airfoil . 

Photos  were  taken  using  a  Nikon  35nim  camera  and  Kodak 
TMAX-400  ASA  black  and  white  film.  The  shutter  speed  was  set 
to  1/125  seconds  with  an  aperture  setting  of  4.0  for  the  light 
conditions.  Film  developing  time  was  optimized  at  9  minutes  at 
75  degrees  F. 

E.  RESULTS  AND  DISCUSSION 

The  result  for  the  tunnel  off  condition  flow  visualization 
experiment  was  as  expected.  The  wing  in  fact  accelerated  the 
smoke  in  its  vicinity. 

The  result  of  the  additional  rake  flow  visualization 
experiments  are  shown  in  Figures  4.5-4.14.  Figure  4.5  shows 
the  stationary  airfoil  at  zero  degree  AOA.  The  Reynolds  number 
(based  on  airfoil  chord)  is  10,000.  It  can  be  seen  that  the 
airfoil  produces  a  small  wake  with  the  boundary  layer  mostly 
attached.  Figures  4.6  through  4.14  show  the  vortical  wake 
flow  patterns  produced  by  plunge  oscillations  at  various 
frequencies  as  indicated.  Most  of  these  pictures  reveal  the 
propulsive  vortical  street  pattern  discovered  in  Reference  7. 
Previous  experiments  by  Neace,  [ref. 9]  found  that  the  tunnel 
was  too  small  for  the  airfoil  size  used,  but  the  airfoil  size 
for  the  present  experiment  seemed  to  be  optimum,  as  seen  by 
the  long  trail  of  vortices.  The  vortical  patterns  show  that 
the  bottom  vortex  is  rotating  clockwise,  and  the  top  vortex  is 
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Figure  4.6  =  1.8008,  h/2b=,1977,  1.56  ft/s 
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Figure  4.7  =  2.467,  h/2b=. 10204,  1.56  ft/s 
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6.753,  h/2b=.1084,  1.56  ft/s 
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Figure  4.13  Kp  »8.223,  h/2b«. 01913,  1.56  ft/a 
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rotating  counterclockwise,  which  is  a  thrust  producing 
vortical  sheet .  It  can  be  seen  in  the  pictures  that  the 
frequency  greatly  affects  the  vortical  strength  (size) . 
Increasing  the  frequency  leads  to  an  increase  in  wake 
vorticity . 


V.  LIFT  ENHANCEMENT  PRODUCED  BY  A  PLUNGING  AIRFOIL 

A.  THEORY 

Chapter  IV  demonstrated  the  propulsive  capability  of  a 
plunging  airfoil.  The  production  of  thrust  implies  the 
generation  of  a  jet  flow  which,  in  turn,  may  be  utilized  as  a 
boundary  layer  control  device.  Therefore,  an  additional  test 
was  conducted  in  the  NPS  smoke  tunnel  in  order  to  explore  the 
feasibility  of  this  concept. 

B.  SETUP 

The  same  NACA0007  plunging  wing  was  used  as  in  chapter  IV 
with  a  different  driving  mechanism.  The  wing  was  mounted  to 
an  ELECTRO- SEIS  Model  113  Shaker  Table  by  APS  Dynamics,  Inc. 
The  shaker  was  located  below  the  test  section  of  the  NPS  smoke 
tunnel  (Figure  5.1).  The  plunging  airfoil  was  mounted  to 
struts  at  both  ends  to  prevent  excessive  bending  while 
plunging.  The  large  airfoil  is  a  cambered  profile  taken  from 
the  rotor  of  a  full  size  helicopter  (13"  chord,  2”  Thickness, 
modified  NACA  airfoil,  Reynolds  Number  of  52,000} .  The  wing  is 
suspended  from  the  tunnel  ceiling  as  shown  in  the  flow 
pictures.  The  design  allowed  for  full  movement  of  the  big  wing 
to  position  it  in  the  vicinity  of  the  plunging  airfoil. 
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C.  WIND  TD21NEL 


This  study  used  the  Naval  Postgraduate  School's  flow 
visualization  wind  tunnel »  The  tunnel  is  an  open- circuit  one, 
with  air  entering  an  inlet  that  measures  4.5  m  X  4.5  m 
(15'X15').  As  the  air  enters  the  tunnel,  it  passes  through  a 
7. 5 -cm  long  honeycomb-  A  9:1  ratio  square  contraction  cone 
directs  the  flow  into  a  test  section  that  is  1.5  m  X  1.5  m 
(5'X5'),  and  6.7  m  (22')  long,  as  seen  in  Figure  5.2.  The 
flow  is  then  exhausted  into  the  atmosphere  through  a  fan, 
which  uses  variable  pitch  blades  to  control  the  speed  of  the 
flow.  The  speed  control  toggle  switch  is  located  right  below 
the  red  and  green  on/off  switch  located  in  the  left  side  of 
the  tunnel  control  room.  The  tunnel  speed  was  determined  using 
a  digital  manometer  which  was  verified  for  accuracy  (Figure 
5.3)  . 

An  observation  booth  is  located  on  the  side  of  the  tunnel. 
A  glass  window,  1.6  m  X  1.1  m  (5.2'  X  3.4'),  provides  the 
primary  viewing  area  from  the  observation  room  and  a  second 
one,  0.4  m  X  1.23  m  (1.33'X4'),  is  located  in  the  tunnel's 
roof.  The  main  viewing  window  had  sufficient  area  for  most  of 
the  photography,  with  the  top  window  used  for  ill\imination.  A 
circular  tumtaUale  was  located  on  the  floor  of  the  test 
section  [ref. 11]  which  allowed  for  easy  access  to  the  shaker 
table.  The  walls  and  floor  of  the  test  section  were  flat  black 
for  low  light  reflectivity. 
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Figure  5.2  Tunnel  Layout 
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0.662  1.5  57.09  53.92  5.88% 


D.  SMOKE  GENERATION 

The  smoke  was  generated  in  the  Rosco  smoke/fog  machine. 
Many  different  smoke  injection  techniques  were  tried  but  with 
less  than  satisfactory  results.  Smoke  rakes  were  first  tried 
outside  the  tunnel  with  the  tube  number  varying  from  2  to  30 
tubes.  The  tubes  were  inserted  in  the  honeycomb  and  also 
separated  different  distances  from  the  inlet  of  the  tunnel. 
The  rake  was  also  tried  inside  the  test  section  with  very  bad 
results  (smoke  dispersed  immediately) .  Problems  ranged  from 
lack  of  smoke  aind  turbulence  when  enough  smoke  was  present. 
The  Rosco  machine  at  its  lowest  setting  was  producing  a  very 
high  smoke  volume  and  whenever  the  flow  was  restricted  by  a 
smoke  rake  the  smoke  production  went  way  down.  The  final 
technique  adopted  was  very  simple.  The  smoke  output  was  sent 
directly  from  the  machine  to  a  1"  nozzle  which  was  manually 
waved  at  the  entrance  of  the  tunnel  to  make  a  steady  cloud. 
The  steady  cloud  was  gradually  pulled  into  the  tunnel, 
producing  a  thick  smoke  sheet  in  the  test  section. 

E.  PHOTOGRAPHY 

Photos  were  taken  using  a  Nikon  35mm  camera  and  Kodak 
TMAX-400  ASA  black  and  white  film.  The  film  speed  was  set  to 
1/250  seconds  with  an  aperture  setting  of  4.0  for  the  light 
conditions.  Film  developing  time  was  optimized  at  9  minutes  at 
75  degrees  F. 
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F.  EXPERIMENTAL  PR0CEDX7RSS 

The  first  step  involved  a  look  at  the  large  airfoil  to 
verify  normal  flow  patterns  (Figure  5.4  and  5.5)  and  find  the 
AOA  for  initial  trailing  edge  separation  (Figure  5.6) .  Next, 
the  plunging  airfoil  was  placed  in  the  tunnel  by  itself  and  a 
run  was  made  to  verify  the  propulsive  capability  in  the  larger 
tunnel  at  higher  speeds.  As  seen  in  Figure  5.7,  the  airfoil 
produced  a  drag  vortical  flow.  Figure  5.8  and  5.9  shows  the 
propulsive  pattern  of  the  propulsive  airfoil.  Finally,  the  two 
airfoils  were  placed  in  close  relative  position  to  see  the 
interference  effect  between  the  two  airfoils. 

6.  RESULTS  AND  DISCUSSION 

Several  airfoil  position  combinations  were  studied,  as 
shown  in  Figures  5.10  through  5.15.  Figure  5.10  and  5.11  show 
the  plunging  airfoil  located  at  approximately  .65  chord  of  the 
large  airfoil.  Figure  5.12  through  5.15  show  the  plunging 
airfoil  located  at  approximately  .75  chord. 

The  differences  between  the  plunging  on  and  off  condition 
were  not  easy  to  see  with  the  eye  but  pictures  indeed  showed 
some  differences  between  the  two  conditions.  A  shortcoming  of 
this  experiment  was  the  inability  of  the  plunging  airfoil  to 
run  parallel  with  the  large  airfoil.  Additionally,  sizing  of 
and  relative  positioning  of  the  two  airfoils  was  not  optimized 
to  give  best  results.  The  two  airfoils  were  chosen  from  the 
resources  available  and  time  constraints  prevented  a  more 
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airfoil  at  zero  AOA 
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Figure  5.10  Large  airfoil  at  12  degrees  AOA,  steady  plunging 
airfoil  position  1 
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Figure  5.11  Large  airfoil  at  12  degrees  AOA,  plxinging 
airfoil  Kpa3.42,  position  1 
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Figure  5.12  Large  airfoil  at  12  degrees  AOA,  plunging 
airfoil  steady,  position  2 


Figure  5.13  Large  airfoil  at  12  degrees  AOA,  plunging 
airfoil  Kp»3.42,  position  2 


Figure  5.14  Large  airfoil  at  14  degrees  AOA,  pliinging 
airfoil  steady,  position  2 
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Figure  5.15  Large  airfoil  at  14  degrees  AOA,  plunging 
airfoil  Kp=3.42,  position  2 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  SINGLE  AIRFOIL  ANALYSIS 

The  modified  version  of  U2DIIF  (UPOT)  can  perform 
aerodynamic  calculations  over  any  range  of  reduced 
frequencies-  The  nonlinear  theory  presented  here  for  harmonic 
motion,  and  the  phase  relationships  that  exist  between  the 
airfoil  motion  and  the  aerodynamic  forces  have  been 
extensively  verified  by  comparison  with  Theodorsen's  linear 
theory.  Furthermore,  this  panel  code  was  applied  to  the 
analysis  of  incompressible  bending- torsion  airfoil  flutter. 
Again,  excellent  agreement  with  the  classical  Theodorsen 
analysis  was  obtained. 

Access  to  faster  con^utational  means  is  recommended  to 
shorten  the  time  needed  to  predict  the  flutter  points.  The 
code  should  be  modified  to  incorporate  three-dimensional 
calculations  which  would  help  solve  more  difficult  flutter 
problems . 

B.  FLOW  VISUALIZATION  EXPERIMENTS 

The  flow  visualization  experiment  successfully  showed  the 
development  of  thrust  produced  by  a  plunging  airfoil.  The 
enhanced  lift  experiment,  on  the  other  hand,  was  not  a 
complete  success.  The  smoke  visualization  presented 
difficulties  that  were  not  satisfactorily  overcome.  As  a 
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result,  the  pictures  taken  were  somewhat  inconclusive. 
Futhermore,  the  angle  of  attack  of  the  oscillating  airfoil 
could  not  be  changed  thus  making  it  difficult  to  achieve  a 
flow  condition  conducive  to  lift  enhancement. 

It  is  recommended  that  further  experiments  be  conducted  in 
the  low  speed  smoke  tunnel  with  a  shaker  table  capable  of 
moving  an  airfoil  at  harmonic  frequencies  near  40  HZ. 
Additionally,  the  airfoil  must  be  modified  to  allow  change  of 
AOA.  Finally,  the  Rosco  smoke  machine  output  volume  must  be 
modified  to  permit  much  lower  smoke  output.  This  final  point 
proved  to  be  the  single  largest  detriment  to  the  visualization 
experiment . 
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